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Abstract 



^^. This article provides a unified treatment of an extensive category of non- 

r~| ! linear classical field models whereby the universe is represented (perhaps as a 

brane in a higher dimensional background) in terms of a structure of a mathe- 
matically convenient type describable as hyperelastic, for which a complete set 
of equations of motion is provided just by the energy- momentum conservation 
C^ I law. Particular cases include those of a perfect fluid in quintessential back- 

grounds of various kinds, as well as models of the elastic solid kind that has 
been proposed to account for cosmic acceleration. It is shown how an appro- 
priately generalised Hadamard operator can be used to construct a symplectic 
structure that controles the evolution of small perturbations, and that provides 
a characteristic equation governing the propagation of weak discontinuities of 
diverse (extrinsic and extrinsic) kinds. The special case of a poly-essential 
model - the k-essential analogue of an ordinary polytropic fluid - is examined 
and shown to be well behaved (like the fluid) only if the pressure to density 
ratio w is positive. 



1. Introduction. 

As a generalisation of the category of media that are elastic in the ordinary varia- 
tional sense 12] , the extensive category of models referred to here as "hyperelastic" 
is characterised by an action density L that can be formulated as a non-linear function 
just of a set of scalar p + 1 scalar fields (y9° , yj^ , ...(^^ and their gradient components 
Lp^^.Lp'^^, ...LpP^ with respect to coordinates x" (a = 0,l,...p) on a p -|- 1 dimen- 
sional worldsheet with codimension g > in a background space time endowed with 
Lorentz signature metric that has components g^^ with respect to coordinates x^, 
(/i = 0,l,...,p + g). 

This work is intended for the treatment of scenarios of the usual cosmological type 
with space dimension p = 3, and it is set up (using a background tensor formalism 
of the kind [Sj that was originally developed for the treatment of conducting cosmic 
strings jH]) in such a way as to be applicable not just to models of the traditional 
kind in which the background spacetime dimension is 4, so that the codimension q 
vanishes, but also to models of more exotic "brane world" varieties, in which the 
background is of higher dimension, 5 or more. 

In a typical cosmological application of such a model, ^9° would represent a 
"quintessence scalar" of the kind commonly invoked [H], to account for the appar- 
ent observation of cosmic acceleration while the other scalars y^^ , (/^^ , ^p^ would be 
interpretable as comoving (Lagrange type) coordinates of a material medium of the 
"normal" kind which in the simplest case would be of ordinary perfect fluid type. 
However instead of being a perfect fluid, the "normal" matter characterised by such 
comoving coordinates could just as well be an elastic solid of the kind envisaged by 
Bucher and Spergel [3 El El- 

In a non-cosmological application for which a model of this hyperelastic kind might 
be used, the "normal" constituent could be that of a solid neutron star crust, within 
which a freely flowing superfluid neutron current would be characterised by the scalar 
field, whose gradient would be the neutron momentum covector, Ha = (p°a- ^"^ that 
case the scalar would have to be of ignorable type (meaning that the Lagrangian would 
depend just on its gradient but not on its undifferentiated value (f°) so that, according 
to equation PHj) below, the neutron current would automatically be conserved. (It is 
however to be remarked that in a realistic treatment of a neutron star crust it may 
be necessary to allow for the possibility ^U] that the superfluid neutron constituent 
may not be separately conserved, so that a more elaborate kind of model would then 
be needed.) 

Whatever its dimension, the background spacetime metric will induce a corre- 
sponding worldsheet metric with components 

9 ah = 9 tJLU^ ^a'^ fi ^ I-'-/ 

and with determinant \'g\ in terms of which the action integral will be expressible as 
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f L\\gf/^(F+^x. (2) 



In order for the system to be considered as regular, the induced metric on the 
worldsheet must of course itself have a Lorentz signature, and furthermore the scalar 
field gradients must be linearly independent so that the fields themselves will be 
adoptable as an admissible set of worldsheet coordinates, for which we shall simply 
have x°' = ^9°, which entails that the Lagrangian density L will depend just on the 
undifferentiated position coordinates ip°' and on the corresponding set oi p{p + l)/2 
induced metric components gab- 

A subcategory of particular interest consists of models for which the Lagrangian 
contributions from "scalar" and "normal" parts separate as a sum of the form L = 
Lg + L^ , in which the "normal" part L-^ is independent oi (p° and of ^^q , while the 
"scalar" part is a (non-linear) function only oi ip° and of its squared gradient as given 
by the single induced metric component 'g^g. 

2. Dynamics of a regular hyperelastic system. 

Whichever coordinate system is used, the worldsheet stress energy density tensor 
will have components given jH] by 

ogab 

and there will be a corresponding world hyper-elasticity tensor defined on the world- 
sheet in the manner introduced by Friedman and Schutz |2j as 

^abcd _ ||-||-l/2 ^Vll^ll -^ J _ rfcdab ( a\ 

ogcd 

which (modulo a proportionality factor of -2) is interpretable as a relativistic ex- 
tension of an ordinary (purely spacelike) Cauchy type elasticity E"'^'^'^ such as will be 
discussed below. 

These worldsheet tensors will map naturally into corresponding background space 
time tensors given by the expressions 

BT 

T"" = T'^'^x'^Xb = 2 7^ + Lg'"', (5) 

(fl^vpa ^ (rabcd^nu p a ^ _^ ^ IT'^^^''" , (6) 

' ' ' ' OgpcT ^ 

in which ^^'^ denotes the (first) fundamental tensor of the worldsheet as defined |S] 
in terms of the contravariant version 'g"'^ of the induced metric by the formula 

r" = r'x>',a^:b ■ iv 

When the codimension q vanishes the overline is redundant here as (I7|) will then just 
give back the contravariant version g'^'^ of the background spacetime metric, and the 
corresponding mixed version ^'^ will then be just the same as the Kronecker unit 
tensor 5(^, but in a background of higher dimension n > p +1 this mixed version g'^ 



will be a non trivial (rank p +1) projector mapping vectors onto their world sheet 
tangential parts, and giving a corresponding world sheet gradient operator 

V. = g^lV, , (8) 

in which the distinguishing overline is again redundent in, but only in, the case of 
vanishing co dimension q. 

As in more general brane models j3] , when the variational field equations ensuring 
invariance of the action with respect to localised perturbations of the worldsheet and 
the dynamical fields theron are satisfied, it will automatically follow as a Noether 
identity that the stress energy given by Q will satisfy a divergence condition of the 
standard form 

v^n = o. (9) 

What distinguishes hyperelastic models from others of a more general kind is that 
in the hyperelastic case no other evolution equations are needed: by itself (0) is not 
only necessary but will also be sufficent to ensure that the variational field equations 
are all satisfied. 

The way this works is that, to start with, the evolution of the world sheet location 
(which will only be needed if the codimension q is non-zero) will be governed as always 
[3] by the orthogonal projection of © which will take the standard form 

T>^''K^f = 0, (10) 

in which K^if is the second fundamental tensor as defined [3] by 

K,/ = ru%9'a- (11) 

As well as the Weingarten integrability condition Kfj^^'' = K^,^'' this defining relation 
ensures the worldsheet orthogonality condition K^y'^'g" = 0, which evidently entails 
that -fC^i/^ itself will vanish in the traditional g = case, for which there are no 
external dimensions so that g"'^ will just be the identity matrix 5^. 

When the condition (jlUp (involving q independent component equations) has been 
satisfied, the remaining part of Q will consist just of its tangentially projected part, 
namely the set oi p +1 internal component equations given by 

r,%T^. = , (12) 

which is just what is needed to determine the evolution of the p + 1 independent 
worldsheet scalar s (f"'. 

3. Symplectic perturbation currents and characteristic equation. 

The sufficiency of the divergence condition Q as a complete set of equations of 
motion can be understood as a consequence of the feature that with respect to a 
reference system of the prefered kind in which the scalars </>" are used directly as 



worldsheet coordinates, the configuration of tlie system will be fully determined just 
by the specification of the background coordinates x'^ as fonctions of these scalars. 

It follows that, with respect to such a preferred system, a perturbation of the 
configuration will be fully determined just by the specification of the corresponding 
background coordinate displacement, 6x'^ = ^'^ say, which, to be dynamically admis- 
sible, must of course be such as to satisfy the linear evolution equation obtained as 
the first order perturbation of ©. When another particular solution, Sx'^ = r/^ say, 
is already available (for example as a trivial perturbation generated by a symmetry 
of the system) then the linear equation governing a generic perturbation ^'^ will be 
conveniently expressible as the conservation, 

V^O^ = , (13) 

of a symplectic worldsheet surface current (such as is also of interest jTTj for the pur- 
pose of quantisation) of the kind that has recently been shown to be straightforwardly 
constructable for a widely extended category of brane systems jl]. As in the case of 
applications |3] to conducting cosmic strings, so also in the more general hyperelastic 
systems considered here, the possibility of expressing the perturbation just in terms 
of a displacement vector ^'^ allows the symplectic current fi" to be given exphcitly by 
an expression of the form 

Q-'il f]} = r]^£)/{^} - eO/{r?} , (14) 

in terms of what I shall refer to as the hyper-Hadamard operator, which is a linear 
differential operator whose action on the vector field ^ is given by a prescription of 
the form 

^/{e}=Vp'^V.e- (15) 

in terms of a corresponding hyper-Hadamard tensor that can be seen j3] to be given 
in terms of the hyper Cauchy tensor (jHI) by the formula 

^,V = 9,pT'''' + 2<t,y. (16) 



It is to be remarked that the standard decomposition 

Qfip = QtMu + -L^i. , (17) 

of the background metric into respectively worldsheet tangential and worldsheet 
orthogonal parts 'g^i, and _L^i, (of which the latter will vanish when the codimension 
q is zero) will engender a corresponding decomposition 

Vp'^=Vp^+^Mr' (18) 

in which the world sheet tangential part is given by 



and the remainder in ()18|) is given simply by 

f,;^7 = ±,,r-. (20) 



As in the simple elastic case ^2], one can obtain the characteristic equation govern- 
ing the propagation of a discontinuity across a worldsheet hypersurface with normal 
covector 

A„ = A^x';,, A^±t = 0, (21) 

of the second derivative of the perturbation vector ^, using the Hadamard rule to 
the effect that it must be given by an expression the form 

[V^V.e] =A^A,C^ (22) 

in which the vector (^ is a measure (whose calibration depends on the normalisation of 
the - possibly null - characteristic covector Aa) of the amplitude of the discontinuity. 
It can be seen to follow from this rule that the discontinuity of the divergence of the 
symplectic current (fT^ will have the form 

V^n^il v}] = v'^.VKKC ■ (23) 

Since the the conservation law (jlHj) is applicable to an arbitrary reference pertur- 
bation solution ?7 (for which undifferentiated components i]'^ may be chosen without 
restriction at any single given point) it can be seen that the characteristic covector 
A^ must satisfy the condition 

^,%''XAX' = 0- (24) 



For extrinsic perturbations of the worldsheet location, as obtained by taking (^ to 
be worldsheet orthogonal, it can be seen from (j^Uj) that we simply recover the result 
(which is already known j2] to apply for a generic, not just hyperelastic, brane model) 
that the characteristic covector must be a null eigenvector of the stress energy tensor: 

g,uC = ^ T^%A, = 0. (25) 

On the other hand for tangential perturbations - the only kind that can exist when 
the codimension q is zero - the characteristic equation will be expressible in purely 
worldsheet tensorial form: 

e=X^,,C => QacC = 0, (26) 

with 

Qac = Qca = ^a'>^bXd ■ (27) 

The condition for Aq to be an intrinsic characteristic covector is thus that the deter- 
minant \Q\ of this symmetric matrix ()27|) should vanish. 
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4. Standard flow decomposition. 

To qualify as hyperelastic in tlie strictest sense, the system should include a sub- 
system of ordinary elastic type, as characterised by the requirement that all the scalar 
fields except one, v^o say, should have spacelike gradients and therefore timelike world- 
sheets that intersect on a congruence of timelike worldlines, whose unit world tangent 
vector, with components u"" say, will be specified, for a = 1, ...,p, by 

nVfa = 0, ^,,nV = -l. (28) 

Subject to the understanding that the small letters refer to an arbitrary worldsheet 
coordinate system, but that the capitals refer to a system of the preferred type spec- 
ified by setting x° = y?" and x^ = V'^, the induced metric components on which the 
action depends can then be listed in a system of the latter type as 

r =r'v\av % , r ° = r'^^v % , (29) 

and 

gAB = g-b^A^^B^ = 7«Vfa^"a = ^^ , (30) 

in which the space projected part of the metric is specified in of the usual way as 

^ab ^ -ab ^ ^a^b ^ ^^^^ 

The generic action variation will be expressible in terms of these quantities in the 

form 

dL dL 

SL = ^6^'> + ^6ip-+ L,, 5-r' + 2L,„ 8^° + L,^ 6^^ , (32) 

which provides the components with respect to the preferred system of the worldsheet 
tensor 

whose contravariant version L"'^ = g°''^g^'^Lcd provides the partial derivatives 

9L .ab 



-L"" (34) 



Jab 



that are needed for the evaluation of the expression Q for the stress energy tensor, 
which works out as 

At the next order one obtains the tensor 



-'abed 



^cd QfgabQfgcd 

which provides the partial derivatives 



(36) 



dgab&9cd 
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j^abcd ^ j^a{c-d)b ^ ga{cj^d)b ^37^ 



that are needed for the evaluation of the expression (JH) for the hyper-Cauchy tensor, 
which is thereby obtained as 

According to p9|) the hyper-Hadamard tensor will therefore be given (using square 
brackets for index antisymmetrisation) by the expression 

j^ahcd ^ ^^abcd ^ 2L'^y^ + 2L'^[V^^ + 'ig^^'^L^^" + L^'V'" , (39) 

from which it can be seen that the characteristic matrix (jTFjl will be obtained as 

Qac = 2{Lacf'' + 2L^\^)hXd. (40) 

The expansion (j32|) provides a corresponding Eulerian variation - as carried out 
for an undisplaced worldsheet location at a fixed value of the background coordinates 
x'^ and metric g^j^y - that will be expressible in the form 

5^L = ^ ^^°+ ^ ^v'+ j:^v%+ Ji^K ' (41) 

in terms of current vectors given by 

which, according to the variational field equations, will have to satisfy worldsheet 
divergence conditions of the form 



5. The separated case. 

In many of the applications of interest a substantial simplification will be provided 
by the separation of the action density as a sum of the form 

L = L3 + L, (44) 

in which the "normal" part, L^^ is independent of (y9°, and (/9°^ so that it has the 
form of the action density of an ordinary elastic medium as specified [H [T2] as a 
function just of '^^ and ^^^, while the remainder Lg is just the action of a simple, 
albeit non-linear, scalar field model, means that it depends only on </?" and on the 
magnitude, /i say, of its gradient 1-form 

/"a = ^%, (45) 



as given by 

F° = r'f^af-ib = -/i' . (46) 

In such a system, neither part will have any dependence on ^-^^ , so the cross terms 
(providing the effect known as "entrainment" ) in the currents (j42p will vanish, 

^.0 = 0. (47) 

The two subsystems will thus be effectively decoupled, except to the extent that they 
interact via their effect on the worldsheet location itself whenever the codimension q 
is non zero, or alternatively, if the codimension is zero, is zero, via their gravitational 
coupling, which can in that case can be be easily incorporated (by requiring the 
background field g^i, to satisfy Einstein's equations or some generalisation thereof) 
but which is not so easy to deal with in a higher dimensional background (for which 
suitable methods of regularisation jTHlin] may be needed). 

The other components of the tensor Lab in the separated system will be given by 

^dL^^ _dL^ . . 

'° dgoo 9(/i2)' ^^""^ 

and by 

L.s = ^,=\{Kl.s-P.s). (49) 

where 7^^ is the covariant inverse of the contravariant base space metric 7-*^ = ^-^^ 
defined by ()3()|1 . and P ^^ will be interpretable as the correspondingly index lowered 
version of the pressure pressure tensor P^^ of the medium, from which its elasticity 
tensor will be obtainable in the usual way P as 

QPAB 

^""c.=2^=^-P^^7o.- (50) 

This means that with respect to arbitrary worldsheet coordinates the pressure tensor 
Pah = P AB^\'^\ °f the elastic medium will have a contravariant version expressible 
as 

p-b = L^^-b_^^ac^M^^^ (51) 

while the ordinary elasticity tensor of the medium will be given by 

^abcd ^ ^^^^ (^^'^b^bc _ 2^«(c^d)fe) + pa{c^d)h ^ ^a{c pd)b _ A-f^^^^'^^^^Y'' Lefgh ■ (52) 



It follows from (j35j) that the total stress energy tensor will be given by the sum 
in which the scalar field contribution will be given by 



T;^ = L3 r" - 2Lf , Lf = -L>V, L'=—^, (54) 



which means that the scalar field will have a pressure Pg and rest frame energy 
density Pg given by 

while the contribution from the "normal" part will be given by an expression of the 
standard form 

T^ = p^u-u'' + P-\ Pn = -^n- (56) 



6. Separated characteristic equations. 

Like the stress tensor, so also the hyperelasticity tensor will be expressible in the 
separated case as a sum 

rf-abcd ff-abcd , ff-abcd (^^\ 

in which, according to ()38j) . the scalar field contribution will be given by 



with 

dV 

while for the contribution from the "normal" elastic medium it can be seen that we 
shall recover the Friedman Schutz j2] formula 

^f^abcd ^ _^Eahcd^pa[c^d)^b^^a^{cpd)b_ il^p^b^^^d^^a^bpcd^^^^^-ab-cd _2gaic^)b^ ^ 

(60) 

In the separated case the crossed element, with respect to the preferred coordinate 
system, of the characteristic matrix (J27jl will automatically vanish 



Qo. = , (61) 

and there will be a decoupling of the "normal" characteristic modes, for which the 
discontinuity amplitude C" is orthogonal to /i^, from the "scalar" characteristic mode 
for which C^ is aligned with u"". The latter is given by 

e = ^ Qoo = 0, (62) 

so that one obtains a scalar characteristic equation in the form 

(2L;'(/i'^Aj2 -L'^X\ X'= r'')KXb = . (63) 

As the relative propagation speed Wg will be definable using a standard calibration 
of the characteristic covector Xa by 

v^ = (P'^A,) V/i^ , A^ = 1 - ^^ (64) 
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it can be seen that the characteristic equation just gives the condition 



v;^ = 1 + 2fi^L';^/L'^ . (65) 



It can thus be seen that for a scalar model of the "standard" kind with a linear 
equation of state, meaning one with L'' = 0, the discontinuities will travel at the 
speed of light, v^ = 1, while for other models of k-essence fHl or more general kinds 
P^ it can be seen that the condition for causality, v^ < 1, and the reality condition 
for local stability, v^ > 0, will both be satisfied if and only if the scalar equation of 
state is such that 

2fi^L''/L'>0. (66) 



S' 



As well as these simple "scalar" characteristic modes (not to mention the extrinsic 
characteristic modes that may exist, subject to ()25|1 . if there is a higher dimensional 
background) there will be different kinds of "normal" characteristic modes given by 

C° = o => Q^cC^ = o, (67) 

so that using the standard calibration 

K = l^a + VUa , UaU" = , U^Ua = 1 , (68) 

to define the relative propagation velocity v and propagation direction i/° (modulo a 
sign that may be chosen to make v positive) the corresponding characteristic condition 
on the normal covector Aq will be the vanishing of the p-dimensional determinant of 
the matrix whose components can be seen from ()40|) to be given by 

Q,^ = 2iL,^\' + 2L,^^^u-u-) , \' = l-v\ (69) 

For the simple fluid case - as given by the restriction that L^^ should depend only 
on the components of the base metric 7^^ only via its determinant I7I - one will 
obtain 

SO the pressure tensor will have the isotropic form, 

pab^p^^ab^ (71) 

in which the ordinary pressure scalar will be given by the well known (but in my 
previous review [12] miscopied) formula 

P, = Z^+2L^. (72) 

In terms of the corresponding bulk modulus, namely 

/?N = -2P^ = -2|7l|^ = -2L^-4L^ (73) 
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one will obtain 

^ABCD lP 1 AbI CD ~r 2 n'^ Acl BD "•" ^'^'^ Ayo"^ B]C ' \'^) 

and the elasticity tensor will be given ^Hj by 

^abcd ^ ^^^ _ p^^^-b^cd ^ 2P^Y^'Y^'' ■ (75) 

It can thus be seen that - as in the more general case of an isotropic solid configu- 
ration P - the characteristic equation ()67|) for the "normal" fiuid will have solutions 
of two kinds, namely longitudinal modes with v = vi^, and transverse modes with 
V = vs, of which the latter are non propagating in the fiuid - as opposed to solid - 
case, 

C^i.^ = ^ v^ = vi, vl = 0, (76) 

while the former are just ordinary sound waves, 

C- = u- => v' = vl, vl = PJ{p^ + P^). (77) 



7. Discussion: polytropic and poly-essential equations of state. 

For both a "normal" fiuid and a "scalar" constituent, an important role is played 
by their respective pressure to density ratios, namely 

wn = ^n/Pn ) ^s = PsIPs ' (78) 

which can be seen to be given by 

WN = l + 2LyL,, w^^ = 2^i^L'jL^-l. (79) 

In recent years it has become common, in specialised cosmological (though not general 
astrophysical) literature, to use a rather loose terminology whereby the pressure to 
density ratio is referred to as the "equation of state" , an appelation that is jusifiable 
only if the ratio in question is actually constant. If - as will often but not always 
be a good approximation - it can be supposed that this ratio, w^ or ws as the case 
may be, is constant then it will indeed characterise a corresponding equation of state. 
In the "normal" fiuid case, such an equation of state will be of what is known (in 
the astrophysical literature) as polytropic type, with polytropic index 7^, as given an 
ansatz of the form 

^N = -CItI"^^/' , 7N = l + w^N, (80) 

for some constant coefficient C . In the "scalar" case, the postulate that ws should 
be constant can be seen to imply that the equation of state will be of what may be 
termed poly-essential type, as given by a power law ansatz, with index as, of the 
form 

Lg = j^/i"^ as = l + Ws\ (81) 
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in which the coefficient ^ is independent of the field gradient magnitude /i, but is 
given as some function just of the scalar field magnitude ip° . This means that the 
poly-esssential ansatz (J8T| characterises a special subcategory within the category 
of quintessential models called k-essential [ini El CH] • This subcategory includes a 
model of the "standard" type, namely the trivial case of a free massless scalar field, 
in the limit when ws = ^i which corresponds to as = 2. 

Attention in cosmology has centered during recent years on the observational 
evidence to the effect that the universe is accelerating, which suggests the need [20] 
for a model with a rather strongly negative value, somewhere in the range between 
-1/3 and -1, for the mean pressure to density ratio w. This poses a problem for a 
"normal" fluid model of the kind specifled for a given value of w^ by the equation 
of state ()80p for which it is well known that the longitudinal modes (f7^ will have 
propagation velocity given by 

vl = wn , (82) 

so that the conditions of reality (for local stability) and causality imply the restric- 
tions 

< m;n < 1 • (83) 

It has been rather unfairly suggested that, compared with such ordinary fluid 
models, k-essential and other scalar models are advantaged by the absence of such a 
restriction, since they have squared characteristic velocity, f | - as given by equation 
(7) of the article fH] referred to - that "can be positive for any" value of the relevant 
"equation of state" ratio, ws- 

The reason why this is unfair is that if ws is a bona flde "equation of state" 
parameter, meaning that it is actually constant, then the corresponding equation of 
state, namely the poly-essential ansatz (jHTj) entails, according to ()65|) . that that the 
relevant propagation velocity will be given by 

vi = {as-l)-' = ws, (84) 

which means that the situation will just the same as in the ordinary fluid case, in so 
much as the corresponding restriction 

< u^s < 1 • (85) 

will have to be satisfled. 

The way suggested by Bucher and Spergel [Zj for getting over the restriction (jH3I) 
is to a seek a mechanism providing elastic rigidity. As an alternative, the way the 
advocates of k-essential (and other quintessential) cosmological models fHl propose 
to get around the restriction (jH3j) is of course to drop the postulate that ws should be 
a genuine "equation of state" parameter, and instead use equations of state of more 
general kinds in which ws is demoted from the status of a flxed parameter to that of 
a variable field. What is unfair is to give the impression that the possibity of doing 
that is a privilege distinguishing scalar field models from "normal" fiuid models: in 
fact the use of more general kinds of ordinary fiuid model - in which the density to 
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pressure ratio w-^ is just a variable field that may be quite different from the squared 
sound speed - is actually commonplace in many areas of astrophysics, and should not 
be prematurely ruled out of consideration in a cosmological context. 

8. Conclusions 

The category of models presented in the preceeding work is sufficient for a wide 
range of cosmological applications of the traditional 3+1 kind in which the codimen- 
sion q is taken to be zero, so that the spacetime metric can be taken to governed by 
the ordinary Einstein equations. 

The treatment here has been set up in such a way as to be applicable also to 
scenarios in which our 4-dimensional spacetime manifold is considered to be a thin 
worldbrane in a higher - meaning 4 + g - dimensional background with a given 
geometry. However for that kind of application it is difficult - despite the efforts of 
many workers in recent years - to see how to allow realistically for the effect of gravity 
except in a Cowling type approximation in which the scale of the perturbations is 
supposed to be sufficiently small (compared with the relevant Jeans length) for their 
self interactions to be neglected. A much studied (albeit only marginally plausible) 
possibility - that is beyond the scope of the present treatment but that does include 
allowance for strong gravitational coupling - is that of the Randall-Sundrum q = 1 
type scenarios ^Tl |22] (and their reflection symmetry violating generalisations [221 
^^) but the unsatisfactory features of such scenarios include the loss of effective 
predictability due to incoming gravitational waves from the bulk. 

In scenarios with codimension g > 2 it is even harder to see how to allow properly 
for gravitation. Attempts to account for the appearance of 4-dimensional gravity 
as a simulated effect j2Sj due to acceleration with respect to a fixed bulk geometry 
tend to predict comportment of scalar - tensor type rather than the pure spin - 2 
gravity that is actually observed. Progress has however been achieved ^J Ej in the 
regularisation of the divergences that will result from true gravity in the bulk, and 
that can be allowed for (if not too strong) within a treatment of the kind presented 
here as an extra contribution to the effective stress energy tensor appearing in the 
extrinsic characteristic equation (J2H)). The condition that the extrinsic propagation 
velocities should always be real and compatible with causality will evidently restrict 
the admissible values of the eigenvalues of the net stress-energy tensor in (j25|) . In 
particular, if it is of the isotropic form 

T^^^ = p ((1 + w)u''u'' + w g^"") (86) 

then the extrinsic propagation velocity ve will be given by 

vi = -w. (87) 

It follows (unless the codimension q is zero) that instead of being subject to a 
positivity condition of the familiar kind exemplified by (jHS)) and (|HH|l the net pressure 
to density ratio w will have to satisfy the negativity condition 

-l<w<0. (88) 
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Such pressure negativity is not incompatible with the observational evidence, but 
does exclude the possibility of describing the universe just in terms of scalar fields 
and fluids that are purely of respectively poly-essential and polytropic type. 
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